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ABSTRACT 

We present a result of a blind search for [Cii]158/xm emitters at z ~ 4.5 
using ALMA Cycle 0 archival data. We collected extra-galactic data covering at 
330 — 360GHz (band?) from 8 Cycled projects from which initial results have been 
already published. The total number of fields is 243 and the total on-source exposure 
time is 19.2 hours. We searched for line emitters in continuum-subtracted data cubes 
with spectral resolutions of ^50, 100, 300 and 500kms“^. We could not detect any 
new line emitters above a Q-a significance level. This result provides upper limits to 
the [Cii] luminosity function at z ^ 4.5 over L[cii] 10® — IO^^Lq or star formation 
rate, SFR ^ 10 — 1000 Moyr”^. These limits are at least 2 orders of magnitude larger 
than the [C ii] luminosity functions expected from the z ^ 4 UV luminosity function 
or from numerical simulation. However, this study demonstrates that we would be 
able to better constrain the [C ii] luminosity function and to investigate possible con¬ 
tributions from dusty galaxies to the cosmic star-formation rate density by collecting 
Cycle 1-1-2 archival data as the ALMA Patchy Deep Survey. 

Key words: galaxies: formation - cosmology: observations - cosmology: early uni¬ 
verse 


1 INTRODUCTION 


The [Cii] ^P 3/2 — >'^Pi /2 fine-structure transition at 

1900.5469 GHz (157.74 /rm) is a dominant coolant of the 
inter - stellar medium (ISM ) in galaxies (e.g., CrawfOTd^£_^ 
1 19851 : IColbert et alj Il999l : iPfAsantos et all l2013l . l2014l) . 
The [C ii] line is one of the brightest lines at far-infrared 
(FIR) and is expected to be a tracer of star formation 
rate in local to dis t ant galaxies (e.g., SmailetalJ_ 20 111 : 
iDe Looze et ah I I 2 OIII , 12014: ISargsvan et aDl20ll[20li) . 


Since the emission line can hold up to 1% of the bolo- 
metric luminosity of a galaxy, it is also very bright, al¬ 
lowing ’blind’ redshifts to be measured. Indeed, ALMA 
observations of sub-millimetre galaxies (SMGs) serendipi- 
tously detected ] C III emission from two SMGs at z = 4.4 
llSwinbank et al.ll2012l) . They obtained the first constraint 
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to the ]Gll] luminosity function (I/[cii] > 10®T©) at z = 
4.4, usi ng the original s urvey area with APEX/LABOGA 
(LESS. IWeifi et ^ I 2009 I ) and the [C ii] redshift coverage in 
the 7.5 GHz bandwidth. However, their constraint was only 
lower limit because their targets are continuum-selected ob¬ 
jects and they could miss continuum-faint ]G II] emitters. 


The ]C II] luminosity function may be a useful 
tool to estimate the cosmic star-formation rate den¬ 
sity at z > 4, and the evolution of star-formation 

rate densi ty can provide tests of galaxy form ation mod¬ 
els (e.g., iBehroozi. Wechsler. fc Gonrov l2013l ). However, 
at z > 4, the contribution from dusty galaxies to 

the cosmic s tar formation r ate density is sti l l un¬ 
certain (e.g.. iB urgyell a et al.l [ 2 OI 3 I: iDowell et all l2014l : 


iMadau fc Dickinson 20141 : ISwinbank et al. 2014 ). In order 
to obtain better constraint to the ]C II] luminosity function 
at z > 4, we carried out a blind search for ]G II] emitters at 
z ~ 4.5 using ALMA GycleO archival data. Throughout this 
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Table 1. Summary of published archival ALMA Cycle0 band? data used for the emission-line search 


Project ID (Region“) 

Targets 

# of Fields 

Antennas 

t'' 

“^exp 

(min) 

RMS= 

(mjy) 

Synthesised beam 
(FWHM) 

Bandwidth 

(GHz) 

Ref.'* 

2011.0.00020.S (EA) 

NGC1614 

1 

16 

51 

0.65 - 0.70 

1".3- 1".5 

7.5 

1 

2011.0.00039.S (EU) 

SBS0335-052 

1 

24 

243 

0.14-0.18 

0".5 - 0".7 

7.5 

2 

2011.0.00097.S (NA) 

GOSMOS 

114 

17-24 

279 

0.67-2.3 

0 

1 

b 

8 

3 

2011.0.00101.S (EA) 

GRBs 

2 

17-20 

91 

0.39 - 0.48 

0''.9- 1".6 

8 

4 

2011.0.00108.S (EA) 

NGC1097 

1 

14 - 15 

59 

0.35 - 0.42 

1".2- 1".5 

7.5 

5 

2011.0.00208.S (EU) 

NGC1433 

1 

19 

101 

0.33 - 0.39 

0".4-0''.6 

7.5 

6 

2011.0.00294.S (EU) 

EGDFS 

120= 

12 - 15 

244 

1.1 - 5.4 

b 

1 

CO 

bo 

8 

7 

2011.0.00467.S (EA) 

VV114 

3 

18 - 20 

85 

0.27-0.30 

0".4-0".5 

7.5 

8 


•^The ALMA Regions, EA: East Asia, EU: Europe, NA: North America. 

^The total on-source time (minutes). 

^T he 1-u sensitivity at 300 km s ~^ spectral resolution . ___ 

^l:llma nishi &; Nakanis Si ll201 3ll, 2:lHu nt et ahl 2014h. 3:IScoyill e et al.l ||2014^. 4:IWang . Chen. &: Huan3 ||201^'1, 5: iFathi et al] ||2013|), 
lizumi et al.l ll2013h. 6: ICombes et al.| ll2013lL 7 : Coppin et alj ||201^'1, ISwinbank et al.l ||2012|. |2014^. iK arim et al.l <l20 13h. iHuvnh et alJ 
<l2013h . iHodge et al.l ll2013h . IWan^^ al.l ll2013h . IPecarl^'^lTll2014li . 8: llono et al.l ll2013lLISliwa et ai.| ||2013h . iTamura et al.l ll2014ll 
®The data of two 2 : = 4.4 SMGs are excluded. 
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Figure 1. Distribution of [Cll] velocity width (FWHM) is shown 
as a function of the [C II] luminosity (-L[c ii]) for a sample of [C II]- 
detected galaxies/QSOs at 2 : > 4 in literature. The JCll] velocity 
width ranges from ~ 50 — 700kms“^. The dashed line shows the 
median JCII] velocity width of FWHM ~ 330 km s“^. The dotted 
lines represent the 10 and 90 percentiles. We use this plot to 
motivate the velocity resolution of the four sets of data cubes 
used in our line emitter search. 
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Letter, we adopt a cosmology with Hq = 72 km s ^ Mpc 
= 0.27, and = 0.73. 


2 DATA AND RESULTS 

In ALMA science archive, we searched Cycle 0 Projects cov¬ 
ering 330 — 360 GHz (in band 7), which corresponds to the 
[Cll] redshift range of z = 4.28 — 4.76. We collected 8 extra- 
galactic projects from which initial results have been already 
published. We only used the published data because we were 
able to know the data quality from the papers before we 
downloaded the data from the archive. We summarise the 
properties of the data sets in Table 1. The data contain 


very deep single pointing and shallow multiple pointings. 
The data with 7.5 GHz bandwidth were taken with spec¬ 
tral mode (FDM mode, the spectral resolution is 0.98 MHz 
or ~0.85kms“^) and the data with 8GHz bandwidth were 
taken with continuum mode (TDM mode, the spectral res¬ 
olution is 31.2 MHz or ~27kms“^)Q The total number of 
fields is 243 after excludin g the data sets of the two [G II]- 
detected SMGs at 2 = 4.4 Iswinbank et aI.|l2012l L The total 
on-source exposure time of the data is 19.2 hours. The pri¬ 
mary beam size is FWHM~ 18” and the synthesised beam 
size ranges from 0”.4 —3”.8. We note that, in our procedure, 
we potentially miss spatially extended [Gll] emitter whose 
size is much larger than the synthesised beam size. 

We used the ALMA data reduction package CASA0 
llMcMulIin et al.|l2007l ') for continuum subtraction and imag¬ 
ing. We used the calibrated data product provided by the 
archive, without any additional re-calibration of the data. In 
order to ensure that the calibrations were correctly applied 
to the data products, we used CLEAN to image the phase 
calibrator from each of the datasets. We confirmed that the 
phase calibrator is detected as a spatially un-resolved source 
at the phase centre and the derived flux density is consis¬ 
tent with the flux densities listed in the ALMA calibrator 
database. We used uvcontsub to subtract the continuum 
from the visibilities using line-free parts of the original tar¬ 
gets with the multiple spectral windows. The continuum 
subtraction procedure is necessary to subtract possible side 
lobes of a bright continuum source and to avoid to detect 
positive noise peaks on the bright continuum source as line 
emitters. We used CLEAN to construct data cubes with nat¬ 
ural weighting to maximise the sensitivity. We adopted a 
pixel size of 0”.2, which is x2 smaller than the smallest syn¬ 
thesised beam size of 0”.4 in the data. 

In order to determine the binning size of the data 
cubes, we checked the [G II] velocity range of high -2 galax¬ 
ies/QSOs in literature. Figure 1 show a distribution of the 


^ In TDM mode, the usable bandwidth is also ~ 7.5 GHz after 
excluding the edges of the spectral windows with high noise. 

^ http://casa.nrao.edu/ 
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Figure 2. Flux distribution of the data cubes with four different spectral resolutions (50, 100, 300, and 500 km s”^) within the primary 
beam before primary beam correction. The flux of each channel is normalised by using the rms. The dotted curves show a Gaussian 
function. The blue dashed vertical lines indicate —6-a and 6-cr, which is used as a detection threshold of line emitters. 


[C II] velocity widt h (FWHM) of 27 [C Ilj-dete cted galax- 


ies/QSOs at 2 > 4 

ICox et al. 2011: Wage et al.l 2010l. 

2012; 

Swinbank et al.ll2012l: Gallerani et al.l 20121: lono et al.l 

200e: 

Carilli et al.l 20131; 

De Breuck et al.ll2011. 2014: Neri et al. 

2014 

: iRawle et al.l 

2014 Wane et al.l l2013l: 

Riechers et al. 

201a 

. 12014: Maiolino et al.l 20051. 20091.12012 

: IWalter et al. 

200c 

; Willott. Omont, & BereerorJ 20131: Venemans et al. 

2012 

). The [C II] velocity width ranges from ~ 50 — 


700kms“^ for a [Cll] luminosity range of 2 x 10®L© < 
^[cii] < 2 X 10'^°L©. The median [Cll] velocity width is 
FWHM 300 km s“^. The 90 per cent of the sample are dis¬ 
tributed in the range of 50 — 500 kms“^. Therefore, we made 
data cubes with four different spectral resolutions (50, 100, 
300, and 500kms“^) to search for line emitters in this ve¬ 
locity range. 

Figure 2 shows flux distributions of the data cubes 
within the primary beam before primary beam correction. 
The flux is normalised by using the rms of each binned chan¬ 
nel. The blue dashed vertical lines indicate — 6 -cr and 6-a, 
which is used as a detection threshold of line emitters. The 
overall shapes of the normalised flux distributions are well 
fitted with a Gaussian (The dotted curves). However, the 
negative tails extend to larger than the 5-cr. Therefore we 
set a threshold of 6 -cr for our line detection conservatively. In 
this search, we could not detect any line emitters above the 
6 -cr significance level. We note that we excluded data with 
emission lines from the original targets and channels with 
very high noise, often seen at the edge of spectral windows. 

We tested the reliability of our line emit ter search as 
follows. Based on the data re-calibrated by iHodee et al.l 


ll2013ll . ISwinbank et al.l ll2012ll detected the z = 4.4 [Cii] 
emission lines at 7.0-cr from ALESS65.1 with L[cii] = 
3.2±0.4x lO®I /0 and at 5.3-cr from ALESS 61.1 with L[c n] = 
1.5 ± 0.3 X lO^L©, respectively. In our procedure, based on 
the data without re-calibration, we detected the emission 
line at 6 . 2 -cr from ALESS 65.1 with L[cii] = 2.9±0.5xl0®L© 
and merginally detected the emission line at 3.8-cr from 
ALESS 61.1 with L[cii] = 1-5 ± 0.4 x 10®L©. Although the 
line luminosities measured in the two different procedures 
agree well, the sig nificances in our proce dure are somewhat 
less than those in ISwinbank et ^ J 2 OI 2 I L These differences 
may come from different methods of data calibration and 
reduction. This test demonstrated that our procedure can 
detect [C ii] emitters if they exist in the data cubes. 

The non-detection of new line emitters provides upper 
limits for the [C II] luminosity function at 2 ; 4.5. We listed 

the survey depths, survey volumes, and 1 -cr upper limits 
in Table 3. We calculated the survey volumes as a function 
of survey depth0 We took into account the primary beam 
correction for the calculation because the line sensitivity is 
not homogeneous even in the same channel. We assume that 
line emitters are spatially unresolved in the data cubes. The 
l-cr confidence upper limits on the space densities of [Cll] 
emitte rs are calculated using Poisson statistics by iGehreld 
(Il986ll . We plot the upper limits for the An = 300kms“^ 


® Data cubes with higher spectral resolutions have higher sen¬ 
sitivities for narrower emission lines and therefore have larger 
survey volumes for a given line luminosity limit. 
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Table 2. Constraints for [Cll] luminosity function at ^ ~ 4.5 



Av = 

50 km s ^ 

Av = 

100 km s 3 

Av = 

300 km s 1 

Av = 

500 kms-l 

(lO^Lo) 

ya 

survey 

(cMpc^) 

ATfe 

(cMpc“®) 

ya 

survey 

(cMpc^) 

Nb 

(cMpc-3) 

ya 

* survey 

(cMpc^) 

N>> 

(cMpc-3) 

ya 

survey 

(cMpc^) 

TV*’ 

(cMpc-3) 

0.13 

24 

< 7.8 X 10-2 

12 

< 1.6 X 10-1 

0 

— 

0 

— 

0.25 

97 

< 1.9 X 10-2 

42 

< 4.3 X 10-2 

12 

< 1.5 X 10-1 

6.1 

< 3.0 X 10-1 

0.5 

7.1 X 10^ 

< 2.6 X 10-3 

1.5 X 102 

< 1.2 X 10-2 

57 

< 3.2 X 10-2 

31 

< 6.0 X 10-2 

1 

3.5 X lO^ 

< 5.3 X lO-'^ 

1.4 X 103 

< 1.3 X 10-3 

2.8 X 102 

< 6.7 X 10-3 

1.0 X 102 

< 1.8 X 10-2 

2 

7.7 X 103 

< 2.4 X lO-'^ 

4.9 X 103 

< 3.8 X 10-1 

1.7 X 103 

< 1.1 X 10-3 

7.9 X 102 

< 2.3 X 10-3 

4 

1.1 X lO"' 

< 1.7 X 10-"^ 

8.4 X lO® 

< 2.2 X 10-1 

4.7 X 103 

< 3.9 X 10-1 

3.0 X 10® 

< 6.1 X 10-1 

8 

1.3 X lO"' 

< 1.5 X 10-"^ 

1.1 X lO"' 

< 1.7 X 10-1 

7.4 X 103 

< 2.5 X 10-1 

5.5 X 10® 

< 3.4 X 10-1 


“The survey volumes are calculated after the prima ry beam correction. 
*’1(T upper limits from non-detection llGehrelslll986l l. 
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Figure 3. ALMA constraints for the [Cll] luminosity function at 2 : ~ 4.5. The black arrows show 1-cr upper limits from non-detection of 
line-em itters in data cubes with 300kms~^ spectral resolution. The red arrow indicates a lower limit at 2 : = 4.44 from ISwinbank et al. 
J2^^). The solid blue curve represents the [Cll] lumino sity function expecte d from the 2 : ~ 4 UV luminosity function jYoshid^^^^ 
l2Q06li and the SFR/L[cii] calibration for high- 2 : galaxies ijPe Looze et al.]f2014^ . The dotted curves show th e dil-a uncertaint i es (0.4 dex) 
of the SFR/L^Qiij calibration. We also adopted the dust attenuation in UV, Ajjy = 1.0 at 2 : 4, by iBurgarella et al.l ll2013l) and 

the relationship between UV luminosity density (Luv) star-formati on rate (SFRO bv[Kennicutt &; E vand l|2Q12lL The dashed curve 
is the [C II] luminosity function predicted from a numerical simulation llOkamoto. Shimizu. & ^stiida[l2^^) with the same empirical 
SFR/L^(^jj] calibration. The current upper limits are still at least 2 orders of magnitude larger than the expected [Cll] luminosity 
functions. The green area indicates the survey volume and depth expected from ALMA Cycle 1 and 2 archival data. 
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case as arrows in Figure 3. We note that, in addition to [C ii], 
[O III] 88 /rm, [N II] 122 /rm, [O l] 145 /im, and [N II] 205 /rm at 
high redshift (2^3), and high-J CO at lower redshift (i.e. 
2 ^ 0.4 for J ^ 3) could be observed in ALMA band? 

fe.g.lSwinbank et al. I 2 OI 2 I : ICarilli fc Walterj|2013l : lOno et al] 

[203)- However, the upper limits for the [Cll] luminosity 

function are not affected from the other possible line emit¬ 
ters at different redshifts. 


3 DISCUSSIONS AND SUMMARY 

For comparison, we plot the [C II] luminosity function ex- 
pected from the 2 ^ 4 UV luminosity function from 
lYoshida et al.l ll2006h in Figure 3. We adopted a dust attenu - 
ation in the UV, Auv = 1 at 2 ~ 4 (iBurgarella et alll2013ll . 
the relationship between UV luminosity density (Luv) and 
star-formation rate (SFR), SFRfMrr^ yr~^) = 1-2 x 10^®Luv 
(ergss“^ Hz“^) by iKennicutt fc EvansI ll2012ll . SFR/Ljcn] 
calibration for high -2 galaxies, IokSFR(M(^ y r ~^) = 

—8.52 -I- 1.18xlogL[cn](L0), by iDe Looze et al'l ll2014h . 
We also compared with the [C II] luminosity function 
at 2 = 4 predicted from a numerical simulation 

JOkamoto. Shimizu, fc Yoshidall2014h with the same empir¬ 
ical SFR/L[cii] calibration. The current upper limits are at 
least 2 orders of magnitude larger than the expected [Cll] 
luminosity functions. Given the depth and sensitivity im¬ 
provements with ALMA in Cycle 1 and 2, we expect to be 
able to detect [C II] emitters in blind searches as ALMA 
Patchy Deep Survey (see the green region in Figure 3). 

If we adopt a detection threshold of 5.5 -ct, we could de¬ 
tect 6 emitter candidates. If these sources are real, we would 
have a number density of ]C II] emitters of ~ 2 x 10“^ Mpc“® 
for L[cii] > 5 X 10® Lq. In this case, the [Cll] luminosity 
function (or the SFR density) at 2 > 4 would be ~ 2 orders 
of magnitude higher than those expected from UV obser¬ 
vations. However, the negative tail of the noise distribution 
also extends to -5.5-a and the number of negative sources is 
2 below -5.5-a. If we consider Poisson statistics, the l-cr un¬ 
certainty of the false positive is 2.0^1'^. The number excess 
of the ’’detected” sources is only ~ 1.5-a compared to the 
false positive. Therefore, the 5.5-a detection threshold may 
be too aggressive to detect real source. 

In summary, we carried out a blind search for [C II] emit¬ 
ters at 2 ~ 4.5 using a part of ALMA Cycle 0 archival extra- 
galactic data. In the continuum-subtracted data cubes with 
spectral resolutions of ~50-500kms“^, we could not detect 
any new line emitters above a 6-a significance level. How¬ 
ever, we obtained the first upper limits to the ]C II] lumi¬ 
nosity function at 2 ~ 4.5 over L[cii] ~ 10® — 10^° Lq or 
star formation rate, SFR ~ 10 — 1000 Moyr”^. These lim¬ 
its are at least 2 orders of magnitude larger than the ]Cll] 
luminosity functions expected from UV observations or nu¬ 
merical simulations. Future ALMA Patchy Deep Survey us¬ 
ing archival data will be able to detect [C II] emitters in the 
blind searches and to study the cosmic star-formation rate 
density at 2 ^ 4. 
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